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Abstract 

Background: The role that the peritumoral stroma plays in the growth of tumours is currently poorly understood. 
In this manuscript the morphometric characteristics of basal cell carcinoma subtypes and their associated 
peritumoral stromas are presented. 

Methods: Ninety eight digitized basal cell carcinoma histology slides were categorized as infiltrative, nodular, or 
superficial subtypes, and were analysed using a combination of manual and computer-assisted approaches. The 
morphometric characteristics of the tumour nests and their associated peritumoral stroma were quantified, and the 
presence of a marked immune reaction or elastosis was noted. 

Results: The tumour to stroma ratio was different among each tumour subtype. Elastosis was identified in a 
greater proportion of the infiltrative tumours. 

Conclusions: Quantitative differences exist between the peritumoral stroma of basal cell carcinoma subtypes. 
Future work exploring the relation between these morphometric differences and biochemical variations in 
peritumoral stroma may further our understanding of the biology of carcinoma development. 

Trial Registration: Not applicable. 



Background 

Carcinomas (cancers derived from epithelial cells) are 
accompanied by a specialized peritumoral stroma which 
plays an important but poorly understood role in the 
growth and metastasis of the tumour [1-4]. The peritu- 
moral stroma is composed largely of activated fibro- 
blasts, inflammatory cells and vasculature and engages 
in a complex crosstalk with the developing tumour. This 
interplay appears to be mediated largely by soluble para- 
crine factors [1-17]. 

Modifications to the peritumoral stroma include 
alterations in the secretion of extracellular matrix pro- 
teins and growth factors by fibroblasts. Particularly 
important are the stromally expressed matrix metallo- 
proteinases (MMPs), which may have wide-ranging 
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effects on tumor growth, angiogenesis and metastasis 
[17]. A number of other soluble fibroblast-derived fac- 
tors have been shown to induce carcinoma initiation 
and progression with Wnt-1, HGF and TGF-|3 playing 
particularly prominent roles [18]. Gene expression pro- 
filing of both stroma and tumor in an inducible human 
neoplasia model showed that during tumour growth epi- 
dermal tissue first showed increased expression of genes 
mediating cellular biosynthesis followed by increased 
expression of genes promoting proliferation and lastly 
increased expression of genes involved in extracellular 
matrix remodelling and increased cell motility. In con- 
trast, gene expression by stromal fibroblasts shows a dif- 
ferent pattern with early increased expression of genes 
involved in angiogenesis [13]. 

In addition to the paracrine factors secreted by fibro- 
blasts, inflammatory cells in the peritumoral stroma also 
secrete growth factors, as well as reactive oxygen species 
which may further effect cancer growth [19]. Likewise, 
complex interactions involving activated fibroblasts and 
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tumor-associated macrophages promote angiogenesis 
[10]. 

As a model for the interaction between epithelial tumor 
cells and stroma, basal cell carcinomas (BCC) offer a 
number of unique opportunities to study these processes. 
In addition to being the most common malignant tumor, 
at least in Caucasians [20], BCC are accompanied by a 
well-defined peritumoral stroma. The contribution of 
this peri-tumoral stroma to the growth of BCCs is parti- 
cularly interesting as several studies have found that the 
growth of BCCs is stroma-dependent and limited by host 
immune response. For example, human to mouse mice 
BCC transplants are only successful in T, B, and NK cell 
deficient mice [21-23]. Moreover, the variety of morpho- 
logical patterns of this tumor provides an opportunity to 
study associations between stromal variation and tumor 
architecture. 

Previous work on BCC peritumoral stroma has defined 
some qualitative features of this tissue. Work by Hum- 
phreys et al. has shown an increased number of mast 
cells and dermal dendrocytes expressing CD34 and GPlb 
adjacent to BCC tumor nests, although no differences 
were noted among tumor subtypes [24]. In contrast, 
among the three most common histological patterns of 
BCC growth (nodular, infiltrative and superficial) striking 
differences in the expression of CD44 glycoprotein 
expression molecules have been observed, with infiltra- 
tive tumors showing the greatest expression of CD44 iso- 
forms as revealed by immunohistochemical staining [25]. 
Dingemans et al. suggest that this finding indicates that 
CD44 may play an important role in BCC biology and 
specifically in determining the histological growth pattern 
and degree of tumor invasiveness. This same study 
reported that a -smooth muscle actin positivity (as a mar- 
ker of myofibroblastic differentiation) tended to be great- 
est in infiltrative tumor areas [22]. Thus there appear to 
be qualitative differences in BCC peritumoral stroma 
which are associated with variations in histological 
growth pattern. 

In addition to biochemical factors, some morphometric 
characteristics of BCCs have been explored. Studies 
performed by Dixon et al. [26] and Jacobs et al. [24] have 
investigated the relationship between tumour nest mor- 
phology and tumour aggressiveness. Significant relation- 
ships were found between non-circular nests, smaller 
degrees of peripheral palisading, and tumour aggressive- 
ness. Dixon et al. [26] and Jacobs et al. [27] also found that 
fibrous stroma was associated with less aggressive 
tumours, and that hyalinization was present in the stromas 
of more aggressive tumours. A couple of studies have also 
explored the relationship between nuclear cell morphome- 
try and tumour aggressiveness [28,29]. De Rosa et al. [28] 
claimed that higher nuclear areas and perimeters were 
associated with increased aggressiveness; however Appel 



et al. failed to replicate this finding [29]. Of particular 
note, de Rosa et al. also observed an increased incidence 
of hyalinization in the stroma of the aggressive tumours 
[28]. Considering that previous studies focused primarily 
on the tumour nests, further investigation into the rela- 
tionships between the peritumoral stroma and tumour 
morphometry may be of interest. 

The purpose of this study is to employ advanced 
image analysis techniques to look for patterns in the 
quantitative relationship between different BCC histol- 
ogy patterns and the BCC peritumoral stroma. 

Methods 

Case selection/image acquisition 

Cases were selected from the routine clinical caseload of 
the senior author. Images of 98 hematoxylin and eosin 
stained slides were digitized using a commercial Aperio 
CS-O slide scanner at 80 x magnification, and sections 
containing BCC were stored as JPEG images (1072 x 902 
pixels). To avoid possible confounding effects of biopsy 
site changes and wound healing on the morphological 
assessment of tumour subtype [30,31], we excluded any 
cases that were previously biopsied. We did not include 
any cases with prior biopsies. Do to a lack of recognized 
classification schemes for basal cell carcinoma stroma, 
we did not attempt to sub-classify stroma types. 

The study protocol was reviewed by the Calgary 
Laboratory Services Research Committee and deemed to 
represent a method development project not requiring 
further ethics approval (CLS Study Code RS-11-504). 

Image analysis 

Printed copies of each of the 98 BCC images were initially 
evaluated manually. In each BCC image five neoplastic 
regions of interest (ROI) were identified, and their stromal 
diameters and tumour radii were delineated. ROIs were 
selected based on the following criteria: (i) the area must 
be free of technical artefacts such as variations in staining 
and cutting of slides and must not include borders of the 
tissue sample or image, and (ii) ROIs must include only 
one subtype of tumor (i.e. not transitional or combined 
tumours). When possible each ROI was selected from a 
separate tumour nest. Multiple ROIs were identified in the 
same tumour nest only when the total number of tumour 
nests was fewer than five. The tumour radii and stroma 
diameters were measured directly from the printed image 
for each ROI. A limitation to this approach arose where 
tumour nests were close together and the intervening 
space was entirely composed of peri-tumoural stroma. In 
this instance stroma attributable to each nest was arbitra- 
rily defined as the stroma between them divided by two. 
Finally, the presence of inflammation or elastosis was 
noted, and the tumour subtype for each slide was categor- 
ized as nodular, superficial, or infiltrative. 



Lesack and Naugler BMC Dermatology 2012, 12:1 
http://www.biomedcentral.eom/1 471-5945/1 2/1 



Page 3 of 7 



In addition to the manual evaluation, the open source 
digital imaging and analysis software ImageJ (version 
1.44p) was used to quantify the morphology of the 
tumour nests in each slide. By measuring the circularity 
of the tumour nest associated with each of the ROIs, 
the tumour nest roundness was quantified. The circular- 
ity measurement in ImageJ is calculated as follows: 

4 • tt ■ area 

circularity = — 

perimeter 1 

A circularity measurement of 1 indicates a perfect cir- 
cle, while more elongated shapes result in lower values 
[32]. In order to avoid artificially low circularity mea- 
surements, tumour nests cut off extensively by the 
image edge were excluded from this analysis. 

Visual tumour classification 

Images were classified into tumor subtypes according to 
the definitions provided by Crowson [33]. The 98 images 
were composed of 23 superficial BCCs, 60 nodular BCCs 
and 15 infiltrative BCCs. The presence of both solar elas- 
tosis and an accompanying chronic immune response 
were subjectively scored with absent and mild counted as 
negative and moderate and marked counted as present. 

Statistical analysis 

Data analyses were performed using IBM SPSS Statistics 
vl9. Box plots of the morphometric statistics and bar plots 
of the categorical statistics were created using R v2.13.0. 
Non-parametric tests were employed, as initial calculations 
indicated non-normal distributions. An analysis of var- 
iance was performed using the Kruskal-Wallis test, and 
Mann-Whitney tests were employed when significant dif- 
ferences between the groups were found. Categorical vari- 
ables were analysed using the Chi-squared test. Interaction 
effects between continuous and categorical variables were 
tested using binary logistic regression. 

Results 

Examples of the different basal cell carcinoma subtypes 
included in this study are illustrated in Figure 1. Statisti- 
cally significant differences were found among tumor 
subtypes with regard to median tumour radii (p = 
0.005), stroma diameters (p = 0.008), and tumour to 
stroma ratio (p = 0.000006). Post hoc pairwise compari- 
sons showed different median tumour radii (p = 0.0009) 
and stroma diameters (p = 0.003) between the infiltra- 
tive and nodular subtypes. For each of the tumour sub- 
types different tumour-to-stroma ratios were obtained 
(infiltrative-nodular p = 0.00006, infiltrative-superficial 
p = 0.006, nodular-superficial p = 0.005). There was 
insufficient statistical evidence to conclude that circular- 
ity varied by tumour subtype (p-value = 0.3). The 




morphometric comparisons are summarized in Figure 1 
and Table 1. 

The presence of elastosis was correlated with the 
tumour subtype (p = 0.033), as eight of the fifteen (53%) 
infiltrative, twelve of the sixty (20%) nodular, and six of 
the twenty-three (26%) superficial tumours contained 



Table 1 Morphometric characteristics of basal cell carcinomas 



Variable 


Infiltrative (n = 1 5) 


Nodular (n = 60) 


Superficial (n = 23) 


All (n = 98) 


P1 


P2 


P3 












(Infiltrative-Nodular) 


(Infiltrative-Superficial) 


(Nodular-Superficial) 


Tumour radius 


52.92 (39.56, 


71.83 (55.45, 


56.62 


64.29 


0.0009 


0.06 


0.3 


(microns)-median 


59.22) 


94.95) 


(50.32, 


(53.16, 








(IQR) 






97.92) 


91.24) 








Stroma width 


34.12 (28.44, 


25.72 (20.28, 


30.17 


27.44 


0.003 


0.5 


0.08 


(microns)-median 


40.56) 


30.17) 


(23.24, 


(21.52, 








(IQR) 






44.02) 


34.86) 








Tumour to stroma 


1 .492 (1 .308, 


2.842 (2.023, 


1.894 


2.170 


1.129e-06 


0.006 


0.005 


ratio-median (IQR) 


1.726) 


4.005) 


(1.627, 


(1.692, 














2.691) 


3.330) 








Circularity-median 


0.3808 


0.4475 


0.4580 


0.4266 


0.1 


0.3 


0.97 


(IQR) 


(0.3158, 


(0.3328, 


(0.2616, 


(0.3155, 










0.4112) 


0.5672) 


0.5813) 


0.5652) 









P 1 ] = p-value obtained for the difference between the median variable values between the infiltrative and nodular subtypes. A two tailed Mann-Whitney test was used with a 95% confidence interval 
P2 = p-value obtained for the difference between the median variable values between the infiltrative and superficial subtypes. A two tailed Mann-Whitney test was used with a 95% confidence interval 
P3 = p-value obtained for the difference between the median variable values between the nodular and superficial subtypes. A two tailed Mann-Whitney test was used with a 95% confidence interval 
IQR = interquartile range 
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Nodular 
Tumour subtype 



Superficial 



Nodular Superficial 
Tumour subtype 



Nodular 
Tumour subtype 



Superficial 



Nodular 
Tumour Subtype 



Superficial 



Figure 2 Box and whisker plots for the mean tumour radii (A), stroma diameters (B), tumour to stroma ratios (C), and circularities by 
tumour subtype (D) 



elastosis. The presence of an immune reaction was not 
found to be significant (p = 0.1). The categorical statis- 
tics are summarized in Table 2. Possible interaction 
effects between dermal elastosis (significant in univariate 
analysis) and the continuous variables were explored 
using binary logistic regression. No statistically signifi- 
cant associations were found using this approach (analy- 
sis not shown). 

Discussion 

This study quantified the morphometric relationships 
between basal cell carcinoma subtypes and their 



peritumoral stroma. Manual and computer-assisted 
methods were employed to quantify size and shape sta- 
tistics for infiltrative, nodular, and superficial BCC sub- 
types. The presence of inflammation and elastosis was 
also noted for each subtype. 

Significant differences were found between the median 
tumour radii and stroma diameters of the infiltrative 
and nodular subtypes. The infiltrative subtype was asso- 
ciated with thicker stroma, but smaller tumour radii, 
with a correspondingly smaller tumour-to- stroma ratio. 
Therefore the biochemical differences previously 
described for infiltrative versus other BCC subtypes may 
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Table 2 Presence of immune reactions or elastosis in basal cell carcinomas. 



Characteristic 


Infiltrative 1 


Nodular 1 


Superificial 1 


P1 


P2 


P3 




(n = 15) 


(n = 60) 


(n = 23) 


(Infiltrative-Nodular) 


(Infiltrative-Superficial) 


(Nodular-Superficial) 


Immune 


4 (0.27) 


1 8 (0.30) 


12 (0.52) 


1.0 


0.32 


0.14 


Reaction 














Elastosis 


8 (0.53) 


1 2 (0.20) 


6 (0.26) 


0.0057 


0.17 


0.37 



Chi-squared tests were used except where an individual cell had a value less than five in which case a Fisher's exact test was used 
1 Data are presented as totals and proportions 

PI = p-value obtained for the difference between the median variable values between the infiltrative and nodular subtypes 
P2 = p-value obtained for the difference between the median variable values between the infiltrative and superficial subtypes 
P3 = p-value obtained for the difference between the median variable values between the nodular and superficial subtypes 



be accompanied by quantitative changes in the peri- 
tumoral stroma. 

Another noteworthy finding was that elastosis was 
present in a significantly higher proportion of the infil- 
trative subtype images. One possible explanation for this 
finding is that the presence of extensive sun damage 
contributes to the phenotype of BCC, predisposing to an 
infiltrative growth pattern. Possible mechanisms for such 
a relationship are speculative but may include differ- 
ences in diffusion characteristics or paracrine factor 
expression in sun damaged skin. Similar arguments 
could be made for the presence of a chronic immune 
infiltrate associated with infiltrative BCC: lymphocytes 
could be recruited in response to some characteristic of 
the tumor or they may themselves secrete factors that 
influence the BCC phenotype. 

Unexpectedly, the circularity of the tumour nest was 
not found to influence the tumour subtype. This may 
have resulted from the presence of nodular tumour nests 
within some of the infiltrating subtype slides, resulting in 
higher mean and more variable circularity measurements. 
Alternatively, variance of the circularity measurements 
among tumor nests on the same slide may provide a 
means of distinguishing the tumour subtypes based on 
their shape. Moreover, the use of different shape metrics 
may provide a superior means of distinguishing tumour 
subtypes based on their morphology. These questions 
could be explored in future studies. 

This study was unique in several ways. First of all, a 
greater emphasis was placed on the relationship of the 
tumour morphometry and its stroma. Previous studies 
have focused largely on the tumour itself, therefore 
ignoring an important factor in tumour growth. Another 
benefit was that the shape statistics were quantified using 
an open-source digital imaging and analysis program. 
This facilitates replication, as ImageJ is freely available 
and is capable of providing circularity measurements 
with the default package. A further benefit of using soft- 
ware to quantify the tumour nest shape is that it elimi- 
nates the subjectivity resulting from categorizing the 
tumour shape qualitatively. Further research could 
improve upon this study in a couple of ways. One 



limitation was that rarer BCC subtypes were not 
included. Also, as mentioned earlier, inherent in an 
observational study is the fact that it was not possible to 
determine whether stromal differences influenced the 
tumour architecture or if they resulted from the architec- 
ture itself. Another limitation of this study concerns the 
applicability of our findings to other malignancies. Basal 
cell carcinomas differ from many other malignancies in 
that they exhibit limited metastatic potential and often 
lack features common to other malignancies such as 
aneuploidy [34]. 

Previous studies that have quantified the morphometry 
of BCC have focused on the relationship of the tumour's 
characteristics and its aggressiveness. Furthermore, these 
studies have focused mainly on the tumour, rather than 
its relationship to its stroma. As a result, it is not possible 
to directly compare the findings presented here with 
those of the previous BCC morphometry studies. That 
aside, previous authors have found that aggressive 
tumours have different morphometric characteristics 
than non-aggressive ones. This may be related to the dif- 
ferent morphology found between the different subtypes 
in this study. The results obtained are also consistent 
with the presence of biochemical differences that were 
previously reported. As discussed, Humphreys et al. [24] 
found that the expression of CD44 glycoprotein differed 
amongst the infiltrative, nodular, and superficial sub- 
types, with infiltrative tumours showing the highest 
expression. The expression of CD44 may have contribu- 
ted to the greater stroma diameters and lower tumour 
radii measured for the infiltrative subtypes in our study. 
Humphreys et al. [24] also reported greater a -smooth 
muscle actin positivity in the infiltrative tumours, and 
Papanikolaou [35] demonstrated integrin-linked kinase 
(ILK) overexpression in infiltrative tumours, which may 
also have contributed to the morphological differences. 

Conclusions 

In this study we demonstrated that there are quantita- 
tive differences in the peritumoral stroma among BCC 
subtypes. Future work should focus on integrating these 
morphometric differences with the previously described 
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biochemical differences in stromal microenvironments 
among BCC subtypes. 
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